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Abstract 

Let Q C R^ be a Lipschitz domain and F be a relatively open and non-empty subset of 
its boundary dQ. We show that the solution to the linear first order system 



VC = GC, Clr = (1) 

Ph' 

vanishes if G G L^fi; R (JVxJV)xJV ) and ( G W 1,1 (f2; R N ), which is the case e.g. for square 
integrable solutions £ of |T} and G G L 2 (S1; R' JVxJV - )xJV ). As a consequence, we prove 

III ■ III : C°(^,r;R 3 ) -> [0,oo), ||sym(VuP- 1 )|| L 2 (n) 

to be a norm for P G L°°(fi;R 3x3 ) with CurlP G L p (n ; R 3x3 ), CurlP^ 1 G L 9 (fi;R 3x3 ) for 
some p, q > 1 with 1/p + 1/q = 1 as well as det P > c + > 0. We also give a new and different 
proof for the so called 'infinitesimal rigid displacement lemma' in curvilinear coordinates: 
Let $ G H^fi;* 3 ) satisfy sym(V$ T V*) = for some * G W 1 '°°(fi;R 3 ) n H 2 (f7;R 3 ) with 
QO \ det > c + > 0. Then there exists a constant translation vector a G R 3 and a constant 

00 ■ skew-symmetric matrix A G so (3), such that 4? = + a. 

m 
m 



Key words: Korn's inequality, generalized Korn's first inequality, first order system of partial 
differential equations, uniqueness, infinitesimal rigid displacement lemma, Korn's inequality 
in curvilinear coordinates, unique continuation 

1 Introduction 



X: 

C$ ' Consider the linear first order system of partial differential equations 

VC = GC, Clr = 0. (2) 

Obviously, one solution is £ = 0. But is this solution unique? The answer is not as obvious as it 
may seem; consider for example in dimension N :— 1, G(t) := 1/t in the domain f2 := (0, 1) with 
r := {0} C <9f2. Then ( := id ^ solves However, in the latter example the solution becomes 
unique if G G L 1 (51), which is easily deduced from Gronwall's lemma. Here we can see that we 
will need intcgrability conditions on the coefficient G; for a precise formulation of the result see 
section [5] The uniqueness of the solution to © makes 

|H|| := ||sym(V U F- 1 )|| L 2 (n) (3) 
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a norm on 



C~(fi,r ; ]R 3 ) := {u G C°°(ft;R 3 ) : dist(supp u, T) > 0}, 
C°°(n ; M 3 ) := {u\n : u G C^°(M 3 ;M 3 )} 

for P G L°°(r2;R 3><3 ) with dctP > c+ > if CurlP G U>(f2; R 3x3 ), CurlP- 1 e L«(ft;R 3x3 ) for 
some p, q > 1 and 1/q + 1/p = 1. Here the Curl of a matrix field is defined as the row- wise 
standard curl in R 3 . 

The question whether an expression of the form ([3]) is a norm arises when trying to generalize 
Korn's first inequality to hold for non-constant coefficients, i.e., 

3c>0 VuGH^TjR 3 ) || symCVwP- 1 )!^) > c|u| H i ( n), (4) 

which was first done for ^P" 1 , CurlP G C 1 (H;]R 3><3 ) by Ncff in [7], cf. [T7]. Here Hi(ft,r;R 3 ) 
denotes the closure of C£°(f2,r;R 3 ) in H 1 (il;]R 3 ). The classical Korn's first inequality is obtained 
for P being the identity matrix, see [HI El El 1131 [El E3 ■ The inequality (0]) has been proved in [TTJ 
to hold for continuous P" 1 , whereas it can be violated for P" 1 G L°°(f2) or P _1 G 5*0(3) a.e.. 
The counterexamples, given by Pompe in [TjJ and [18], see also [16], each use the fact that for 
such P an expression of the form of ||| • ||| is not a norm (It has a nontrivial kernel) on the spaces 
of functions considered. Quadratic forms of the type (@| arise in applications to geometrically 
exact models of shells, plates and membranes, in micromorphic and Cosserat type models and in 
plasticity, H EH ED M ■ 

The so called 'infinitesimal rigid displacement lemma in curvilinear coordinates', a version of 
which can be found in [1] and which is important for linear elasticity in curvilinear coordinates 
(see also [UH]) states the following: If f2 C R N is a bounded domain, * g W 1 '°°(ri; R N ) satisfying 
detV* > c+ > a.e. and $ G H 1 (fi;E Ar ) with sym(V$ T V*) = a.e., then on a dense open 
subset of 51 there exist locally constant mappings a : ft — > H N and A : Q — > so (TV) such that 
locally $ = + a. If H is Lipschitz then the terms 'locally' can be dropped. In their proof pQ, 
the authors apply the chain rule to := c&ovl/^ 1 and use the observation that the conditions 
sym(V$ T V v I') = and sym(V<I>(V , I')~ 1 ) = are equivalent by a clever conjugation with (V^) -1 , 
this is 

(V*r T sym(V$ T V*)(V1') _1 = sym(V$(V*) _1 ) = sym(V($ o * -1 )) o (5) 

together with the classical infinitesimal rigid displacement lemma applied on 0, defined on the 
domain ^(0). If to this lemma a boundary condition $ = on a relatively open subset of the 
boundary is added, one obtains $ = (cf. [2j 1.7-3(b)]). 

The main part of our proof for ||| • ||| being a norm is also concerned with obtaining u = 
from sym(VuP _1 ) = 0. By taking P = to be a gradient, we present another proof of 
the infinitesimal rigid displacement lemma in dimension N = 3 which yields $ = + a with 
A G so(N), a G R . We need slightly more regularity but do not use the chain rule for 0. 

The key tool for obtaining our results is Neff's formula for the Curl of the product of two 
matrices, the first of which is skew-symmetric. We state a generalization of this formula in section 

ED 

This paper is organized as follows: The next section states the main results that will be proven 
in the subsequent chapters. Section 13.11 provides a tool that gives £ = on lines and is used in 
section 13.21 where this is extended to cubes. Section 13.31 then takes care of the whole domain if 
a '(£ = 0)-cube' is given as starting point. In section |3~^1 the uniqueness theorem of section 
is proven, mainly by putting together the results of the previous sections. After that and before 
applying the theorem we have a closer look at the formula for the Curl of a product of matrices 
(section |4.1[) . Finally, in sections 14.21 and l5l respectively, we prove that ||| • ||| is a norm and we 
present our new proof of the infinitesimal rigid displacement lemma. 
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2 Results 



Let us first note that by V we denote not only the gradient of a scalar-valued function, but also 
(as an usual gradient row-wise) the derivative or Jacobian of a vector-field. The Curl of a matrix 
is to be taken row-wise as usual curl for vector fields. 

Theorem 2.1 (Unique Continuation). Let f2 C M. N , N eN, be a Lipschitz domain, T be a relatively 
open and non-empty subset of dfl as well as G £ L 1 {fl;R^ NxN ^ N ). If ( £ W 1 ' 1 (fi; R N ) solves 

VC = GC, Clr = 0, 

then C = 0. 

From the differential equation itself it is not a priori clear that £ belongs to W 1,:L (f2). But 
this can be ensured by requiring higher integr ability of G and £, since for bounded domains, 
e.g., the conditions G £ L 2 (fi) and ( £ L 2 (fi) imply VC £ L 1 ^) and hence C £ W 1 ' 1 ^), where 
an application of the theorem ensures f = 0. Thus we have obtained the uniqueness of ^(re- 
solutions if the coefficient G arc squarc-integrable. Of course, the same holds if £ £ L p (f2) for 
arbitrary p > 1. Then G at least needs to be an L 9 (f2)-function, where l/p+ l/<? = 1. 

Theorem 2.2 (Norm). Let SI C R 3 be a Lipschitz domain, ^ T C <9S1 oe relatively open, 
P £ L°°(f2;R 3><3 ) wii/i detP > c+ > 0, CurlP € U>(ft;R 3x3 ), Cui'lP" 1 £ L«(fi;R 3x3 ) for some 
p, q > 1 m£/i l/p+l/q — 1 . TTien 

I • I : C~(fi, T; R 3 ) [0, oo), u H- | sym(V U p- 1 )| L2(f2) (6) 

defines a norm. 

Remark 2.3. In the case of p = q = 2 and for P £ SO(3) a.e., CurlP" 1 £ L 2 (^) is no additional 
condition, since then CurlP e L 2 (fi) ^> CurlP" 1 £ L 2 (fi). fJNTofe tfterf /or P £ SO(3) a.e. 
generally P, CurlP £ L p (fi) is equivalent to P £ W 1 '^), c/. flf J 

Conjecture 2.4. TheoremW^ holds for P £ L°°(f2) with CurlP £ L p (ft) and det P > c+ > for 

some p > 1 or even p > 1. 

Remark 2.5. Since the norms ||| ■ ||| and || ■ |h 1 (0) are rao ^ shown to be equivalent, it is not clear, 

whether the spaces H 1 (f2,r) = Cg°(f2,r) " Hl(n) and Cg°(i7, T)'" ' coincide. However, by \17\[ , these 
norms are equivalent if P £ C (0) wit/i detP > c + > 0. 

Conjecture 2.6. T/ie norms are equivalent if P £ L°°(f2) urci/i CurlP £ L p (f2) and det P > c + > 
/or some p > 1 or even p > 1. 

Theorem 2.7 (Infinitesimal Rigid Displacement Lemma). Let C M. 3 be a Lipschitz domain. 
Moreover, let $ £ W 1,p (0;M 3 ) and * £ W 1 ' 00 (ft;M. 3 ) nW 2 '«(!l;K 3 ) iratt dctV* > c+ > a.e. 
and p,q > 1, 1/p + l/q = 1. // 

sym(V$ T V*) = 

t/ien i/iere exist a £ R 3 and a constant skew- symmetric matrix A £ so(3), such that $ = + a. 

Remark 2.8. When comparing two nearby configurations of an elastic body, namely ^ : Q — > R 3 
and \1/ : Q — > R 3 , following Ciarlet ^ we may always write "J = ^ + $ , where $ : 51 — > R 3 is i/ie 
displacement from VP to TTie respective metric tensors of the two configurations are V^fr V$ 
and VvE^V 1 !/. in terms of the displacement $ io lowest order we have for the ^-linearized change 
of the metric 

[V* T V* - V* T V*]i in , (E , = [(V* + V$) T (W + V$) - V* t V*]h„,$ 

= V$ T V* + W T V$ = 2sym(V$ T V*). 

Therefore, the infinitesimal rigid displacement lemma expresses the fact that if the linearized 
change of the metric is zero, then the displacement must be (the linearized part of) some rigid 
displacement. 
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3 Proof of the uniqueness theorem 

We start with some preliminaries. 



3.1 Vanishing in intervals 

Let — oo < a < b < oo and / := (a, b). 

Lemma 3.1. Let G G L 1 (I;R NxN ), ( G W 1 - 1 ^;^) with (' = G( and ((a) = 0. Then ( = 0. 

Proof. We use Gronwall's inequality. Because C G W 1 ' 1 (/), £ is absolutely continuous and hence 
it can be written as an integral over its derivative: 



C(ar)=C(o)+/ ('(t)dt = ((a)+ G(t)((t)dt => |C(*)| < |C(o)| + / ||G(t)| |C(*)| dt. 
An application of Gronwall's inequality leads to 

|C(x)|<|C(a)|exp( [ X \\G(t)}\dt)=0, x G I, 



which concludes the proof. □ 

3.2 Vanishing in cubes 

Let Q be a cuboid in and let T be a face of Q, i.e., Q = T x I with / from the previous 
section. By pH Th. 2.1.4] we have that for u G L 1 (Q) the following is equivalent: u G \N 1 ' 1 (Q), if 
and only if u has a representative which is absolutely continuous on almost all line segments in Q 
parallel to the coordinate axes and whose (classical a.c.) partial derivatives belong to L 1 (Q). These 
classical partial derivatives coincide with the weak derivatives almost everywhere. In particular, if 
u G \N 1 - 1 (Q) then u 7 := u(j, •) G \N 1 - 1 (I) f.a.a. 7 G T. Of course, the same holds for u G W 1 ' P (Q) 
with p > 1. 

Lemma 3.2. Let G G L^Q; R (Arx7V)xAr ) and C £ W 1 ' 1 ^;^) wtfi VC = G( and (\r = 0. TTien 

C = o. 

Proof. Since 

|C| L i(n)=/ / |C(7^)|dxd 7 = / / |C^(^)| da: d T 

we only have to show C 7 = a.e.. As C G \N 1 ' 1 (Q), C 7 G W 1 ^/) f.a.a. 7 e T by P Th. 2.1.4], 
as mentioned before. Since £' is the last column V^e^ of VC, we have 

C; = VC( 7 , -)e N = G( 7 , 0C(7, -)e N = G( 7) OCr^ == G 7 C 7 . 

For fixed (7, x) G Q, G(7, x) is a linear mapping from S, N to M. NxN , its product with C 7 (x) G M w 
is an element of R NxN and multiplication by e N gives an element of M. N depending linearly on 
C-f(x). Hence, G 1 {x) is a linear mapping from M. N to a.e.. Even G 7 G L (J; M JVxAr ) holds, since 
G G L 1 (Q). Also, C|r = implies C 7 ( a ) = f.a.a. 7 G T. By Lemma \3. II we obtain £ 7 = 0. □ 

3.3 Unique continuation 

Lemma 3.3. Let G G L x (0; r(NxN)xn^ and (■ g W i.i(^ ; M JV) ^ V( - = G £ Moreover, Zei C 
vanish in an open ball Bcfi. TTien C = 0. 
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Proof. Let CI be convex and pick some x\ G B. Then we can take a straight line between x\ and 
some other point xi £ CI and a cuboid <5 containing this line and having one face being entirely 
located in B. Then by Lemma l3.2l (T = in Q and hence in a whole neighborhood of X2- Since X2 
was arbitrary, we have C, = in Cl. By induction this can be carried over to connected unions of 
finitely many convex sets and hence works for path-connected sets, because every path between 
two points can be covered by such a finite union. Since domains arc path-connected, we finally 
achieve ( = in Cl. □ 

Remark 3.4. Lemma \3.3\ can also be stated as: The equation V£ = G£, i.e., the operator V — G 
has the unique continuation property. Moreover, it is enough that £ vanishes on a small part of 
some (N — 1) -dimensional hyper-plane. 



3.4 Proof of Theorem I27H 

Let C be as in Theorem 12. II If we can show that £ vanishes on an open set, we can apply Lemma 
13.31 and hence, £ must vanish in the whole of Cl. To make ( = on an open set, we transform a 
part of r, where we know £ to be zero, and a neighborhood U onto a cuboid Q, where we can use 
Lemma 13.21 and the transformed function is forced to vanish, hence also C must vanish on U. Let 
us pick a point on T and a corresponding open neighborhood U as well as a bijective bi-Lipschitz 
transformation cp : Q := (—1, 1) N — > U, mapping the cuboid Q := (—1, l)^ -1 x (0, 1) onto U D Cl 
and (-M)^ 1 x {0} onto U n L. Now G := G o ip g l 1 (Q;R( NxN ^ xN ) and, see again e.g. [H 
Th. 2.2.2], C := ( o ip e W 1 ' 1 (Q;K JV '). By the chain rule we have 

VC = ((VC) o ip)W(p = ((G() o<p)\7ip = G(\7tp =: G(. 

Since is uniformly bounded we get 

Vz e Q, y e R w ||G(z)y|| < |G(z)|| |^| ||V^(z)|| < c||G(z)|| [ y | = 

and hence 



\G(z)\\<c\\G(z)\\ 



\\G(z)\\dz<c \\G(z)\\dz = c |G(x)| |detV^ _1 (a;)|da; < c / ||G(x)||dx < 00 
Jq J<p{Q) Jnnu 

since dctV^ -1 is uniformly bounded as well. Thus G G L 1 (Q; ^( NxN ) xN y Because C vanishes 
on UDT, C vanishes on F := ^(UnT) = (-1, l)^- 1 x {0}. Hence C £ W^^Q;^) solves 

VC = GC, C|f = 0. 

Lemma 1531 implies C = in Q. Thus £ = £ o i^ -1 = in {/ HO, which contains an open ball. □ 



4 Proof of the norm property 
4.1 Curl of matrix-products 



We identify 



»3x3 



and R 9 by the following isomorphisms: 



mat 



p 3x3 







a 1 


ft 2 


03 




>-> 


n 4 


«F, 


a 6 


_a 9 _ 






ag 


a 9 _ 



vec := mat 



p 3x3 



We also use the following canonical isomorphism to identify R 3 and so (3): 





" 


-«3 


a-2 




ai 


axl : so (3) -> R 3 , 


03 





-01 








_-ft2 


ftl 







. fl 3. 
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Moreover, we define 









"© 


a 2 


a3 




ai 




diagvec : 


R 3x3^ 


R 3 , 


a\ 


© 


rto 




a 5 










_a 7 


as 


©. 




_a 9 _ 










ai 


© 


©" 




-a 6 


skewvec : 


R 3x3^ 


R 3 , 


(7. 4 


0-5 


© 




03 








_a 7 


ag 










R 3x3 _j 




a i 




"3 




a g 




symvec : 


R 3 , 


© 


"r, 


«(> 




-a 7 








_© 


© 


a 9 




(24 



We note skewvec = symvec = axl, diagvec = on so (3). Furthermore, Ax = axl(A) x x and 
axl _1 (a)a; = a x x holds for all A G so(3) and all a, x G R 3 , where x denotes the cross-product. 
For a matrix Y G R 3x3 with Y T = [yi 2/2 2/3] and vectors y n G R 3 we define 



L 



diag, Y 



L 



skcw.y 



sym,y 



- / sym,y 



axl" 1 yi 





axl" 1 2/3 








axl" 1 j/2 

axl" 1 2/3 

axl -1 2/3 axl" 1 y 2 



- axl" 1 2/2 axl" 1 j/i 

- axl" 1 2/3 
axl" 1 j/3 

- axl" 1 2/2 axl" 1 2/1 








axl" 1 2/1 


axl" 1 2/2 
- axl" 1 2/1 




and note Ly = Ly . Furthermore, for vector fields v in R 3 we set 

Vw := vcc Vu, 

denoting the vector-field containing the nine partial derivatives of the three components of v. 

Now, we extend Neff's formula from [7J Lemma 3.7] in two ways, such that it can be applied 
with weaker differentiability and for general matrices. For this, we define 



W s (Curl, ft; 



I>3x3\ 



:= {Y E L s (ft; 



d3x3\ 



CurlF G L s (ft; 



c 3x3 



)}■ 



Remark 4.1. For skew- symmetric matrix fields we have W s (Curl, ft; R 3x3 ) = W 1,s (ft; R 3x3 ) 
since in this case the Curl controls all the derivatives, see F7I 



Lemma 4.2. Let r,s G (1, 00) with l/r + l/s = 1. Moreover, let X G W 1 - r (ft; R 3x3 ) and 
Y G W s (Curl,ft;R 3x3 ). Then XY G W 1 (Curl, ft; R 3x3 ) and 

Curl(Xy) = mat(Ldia g ,Y V diagvec +L s kcw,y V skewvec +I sym ,yV symvec)X + X Curl Y (7) 

with Ldiag.y , iskcw,y , ^sym.y G L s (ft;R 9x9 ). for skew- symmetric X formula ([7]) turns to 

Curl(XY) = mat L Y ( V axl X) + X Curl Y, (8) 
where det Ly = — 2(detY) 3 . Hence, ifYis invertible, so is Ly- 

We note that for smooth (C 1 ) matrices X, Y, where X is skew-symmetric, formula © was 
already shown in pj Lemma 3.7]. 



G 



Proof. Since C°°(f2) is dense in both W 1 ' r (0) and W s (Curl, 0) we have to show {7J only for smooth 
matrix fields. But this is a straight forward calculation, which we present in the appendix. ([8]) is 
a simple consequence from ([7]) and the assertion about the determinants has been proved already 
in [3 Lemma 3.7]. □ 



4.2 Proof of Theorem HH 

Let u G C;f (ft,L;IR 3 ) with |||u||| = 0. We have to show u = 0. Note that sym(VuP _1 ) = implies 

VuP- 1 = A, (9) 

where A is some skew-symmetric matrix field. Moreover, since AP = Vu we have 

Curl(AP) = 0. (10) 

Without loss of generality we assume T to be bounded (otherwise, replace T by a bounded 
open subset of itself) and that the compact set suppu and L are both contained in some open 
ball B. Define tt := fl n B. Then Vu^F" 1 and A belong to L°°(ft,R 3x3 ) C L r (ft,R 3x3 ) for all 
r G [1, oo]. 

Since CurlP" 1 G L«(f2, R 3x3 ), Lemma PI with X := Vu and Y := P" 1 together with Remark 
Oshow A G W 1 ' 1 (fi,M 3x3 ) and by even A G W 1 '^, R 3x3 ) holds. Another application of 
Lemma IO with X := A and Y := P gives by © and §TU$ 

matip(VaxlA) + A Curl P = 0, 

since A is skew-symmetric. Thus, ( := axl A G W 1 ^(il,]R 3 ) C W^^R 3 ) solves 

VC = - mat Lp 1 vec ( axl" 1 C Curl P) =: G P (. (11) 

Since L P ,Lp X G L°°(f2,R 9x9 ) and CurlP G U>(f2,lR 3><3 ) C L 1 (0,M 3x3 ), also G P belongs to 
LP(0,R( 3x3 ) x3 ) C L 1 (f2,R( 3x3 ) x3 ). Additionally, A and hence < vanish on T by © since u does. 
By Theorem 12. 1[ £ and therefore A and Vu vanish in f2. Thus u = in f2 because u vanishes on 
r. Since supp u C O we finally obtain u = in f2. □ 



5 Proof of Theorem 12.7 

A := V$(V V I') _1 G L p (£1; R 3x3 ) is skew-symmetric by <j5j) - Since the standard mollification pre- 
serves skew-symmetry we can pick a sequence (A n ) C C£°(r2; R 3x3 ) of skew-symmetric smooth 
matrices approximating A in L P (H). Applying Lemma i.e., ([5]), to A n V^ we get 

Curl(A„V*) = mati v *(VaxlA n ) 

with invcrtiblc Ly* G W 1 ' <? (f2) satisfying Ly^ G W 1,9 (fi) by assumption on the regularity of 
Pick 6 G C;f (0;R 3x3 ). Then L^vecG G \Nl' q {VL) and since AW = V<J> G LP(0;R 3x3 ) with 
Curl(AV*) = we have 

(A„V*, Curl(matL v ^ vcc6)) L2(n) -> (AV*, Curl(mat L v ^, vec 6)) L2(n) =0. 

On the other hand we have for the left hand side 

(A n W, Curl(mat L v ^ vec 6)) L2(n) = (Curl(A„V*), mat L v l vec 6) L2(n) 

= ^i v *(VaxlA„),Ly^ vcc0 ) L ^ (a) 

= ( VaxlA„, Ly^L^I vec9) 

\ /L 2 (n) 

= (VaxlA n ,6) L2(n) 

= (axlA„,Div6) L2(fi) -> (axlA,Div9) L2(n) . 

Hence, V axl A = and therefore A G so(3) is constant. Thus, V($ - A*) = V$ - AV* = and 
$ = A* + a with some a G R 3 . □ 
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A Appendix 



We show (J7|) for smooth matrix fields X = [x n m]n.m=i,2,3 and Y = [ynm}n.m=i,2,3- The i-th row of 
XY is the transpose of the vector having the entries Xi n y n k for k = 1,2, 3. Thus, the l-th row of 
Cur\(XY) is the transpose of the vector having the entries di(xi n y n j) — dj{xi n y n i) for k = 1,2, 3, 
where the curl of a vector field v is written as 



curl v 



d 2 v 3 - d 3 v 2 
d 3 vi - 9iu 3 
d\v 2 - d 2 vi 



[diVj - djVi\ k = 



1,2,3- 



Therefore, 



[Curl(AF)] ifc = d l xi n y nj - d xi n y ni + xi n (diy nj - d y ni ) 

= [CurlY]„ fc 

= d l xi n y n] - djXi n y ni + [XCurlY]; fc - 
With the transpose of the n-th row of Y denoted by [y n ]j '■= ynj, we get 

[Gxa\{XY)] lh - [XCur\Y] lk = [Wx ln x y n ] k 
and hence for the Z-th row [Cml(XY) — X CurlF]/ = [(Vx;„ x y n ) T ]i- Finally we obtain: 



Curl(XF) -A Curl F 



(Wx u x y x ) 
(Vx 22 x y 2 ) 
(Vx 33 x y 3 ) 

(axP 1 yiVa;ii) T 
(axP 1 y 2 Vx 22 ) T 
.(axP 1 y 3 Vx 33 ) T 



(Vxin x y n ) T 
(Vx 2n x y n ) T 
(Va; 3 „ x y n ) T 

(Vxi 2 x y 2 ) T + (Va;i 3 x j/ 3 ) 
(Vx 2 3 x y 3 ) T 






(V»21 X ?/i) T 

(Va; 3 i x yi) T + (Va; 32 x y 2 ) 



(axl y 2 Va;i2) + (axl y 3 Vxi 3 ) 
i&xr 1 y 3 Vx 23 ) T 








(axl 1 2/iVx 2 i) T 
(axl -1 2/iVa; 3 i) T + (axl -1 y 2 Vx 32 ) T 



— mat 



axl" yi 





axl 1 y 2 






axP 1 y 3 



mat 



mat 



axl y 3 

-axl -1 y 3 




axl" 1/2 -&xr 1 yi 







Va?u 




Vx 2 2 




_Vz 33 _ 




"V 




"-V.X2 








VZ13 








-V:ri 






Va; 32 








Vz 3 i 






VX 2 1 



mat(L c iiag,i'VdiagvecX + L s kew,yV skewvecA + Z/ symj yVsymvec A) 
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